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The Gaussian chain in a quenched random potential (which is characterized by the disorder 
strength A) is investigated in the d - dimensional space by the replicated variational method. The 
general expression for the free energy within so called one - step - replica symmetry breaking (1 
- RSB) scenario has been systematically derived. We have shown that the replica symmetrical 
' (RS) limit of this expression can describe the chain center of mass localization and collapse. The 

critical disorder when the chain becomes localized scales as Ac — b'^ N~^^'^^'^ ( where b is the 
length of the Kuhn segment length and A*' is the chain length) whereas the chain gyration radius 

7?g ~ 6 (b'' / A) ^ * . The freezing of the internal degrees of freedom follows to the 1-RSB - scenario 
and is characterized by the beads localization length C^. It was demonstrated that the solution for 
appears as a metastable state at A = A^ and behaves similarly to the corresponding frozen 
states in heteropolymers or in p - spin random spherical model. 
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PACS numbers: 36.20.-r Macromoleculos and polymer molecules - 05. 40. -a Fluctuation phenomena, random 
processes, noise and Brownian motion - 75.10.Nr Spin-glass and other random models - 71.55.Jv Disordered 
structures; amorphous and glassy solids 



I. INTRODUCTION 



' O The behavior of polymer chains (with or without excluded volume) in a quenched random potential is a simple 

^ ] but nontrivial problem of statistical mechanics with disorder P]. Despite the vast investigations of both statical 
p, 3, 4, 5, 6, 7, 8, 9, Kdi ITU and dynamical 12, 13, 14] properties many aspects remain to be elucidate. First of 
all we emphasize that one should discriminate between the annealed and quenched problems. In the case that the 
effective diffusion coefficient is large enough, so that the chain experiences different disorder realizations, the problem 
is equivalent to the annealed one (even though the random potential is quenched) . In this case the presence of the 
disorder manifests itself (in the statics) through the reduction of the excluded volume parameter v, i.e. v ^ v — A, 
' where A is the second moment of the quenched random potential. Strictly speaking it was shown long time ago 
' mainly by the use of Imry - Ma arguments, that in the infinite volume limit the quenched and annealed problems 
coincide, so that the Gaussian chain collapses to a size of one monomer. Only if the volume of the random medium 
is finite the chain acquires the characteristic gyration radius Rg ~ b[{A/b'') \n{fl/b'^)]~-^/'-^~'^\ where b is the Kuhn 
. segment length. Nevertheless this approach completely ignores the very important dynamic aspects of the problem. 
Indeed, as the disorder strength A and the chain length N increase, the diffusion coefficient D falls down dramatically 
and the chain can not sample the whole random medium over the course of internal degrees of freedom equilibration. 
It was shown 0,0] that in the d-dimensional space the center of mass diffusion coefficient D « exp[— (A/6'*)iV"], 
where is the Rouse diffusion coefficient, b is the Kuhn segment length, a = 2 — fd and v is the Flory exponent. 

■ As a result the time needed to obtain an "averaged" (annealed) statistics grows exponentially with A and N. In a 
P5 ' recent publication |l5j | we have treated this problem dynamically by functional tools based on a Langevin dynamics 
O ■ approach. This approach avoids replicas naturally but by averaging over the quenched disorder inevitably couples 
^ ■ different dynamic trajectories. This leads ultimately to the non - Markovian diffusional slowing down as well as to 
^ the corresponding freezing and non - ergodic regimes for the chain Rouse modes. We have derived the equation of 

motion for the Rouse mode time - dependent correlation function within a Hartree approximation and found that 
^ J this equation has a memory kernel which is generic for the mode - coupling theory of the glass transition 16] . The 

■ self - consistent treatment of the mode - coupling equations leads to the conclusion that the chain center of mass 
" ■ '■ diffusion coefficient decreases according to the law: D sa -Dr[1 — const(A/6'')A^"]. The individual Rouse modes of 

the chain also freeze at a common disorder strength, but their dynamic pathway is strongly selective to the mode 
number. We have recently launched the intensive Monte - Carlo (MC) simulation to support these analytical 
calculations. Indeed we have found that the diffusion coefficient as a function of {A/b'^)N°' drops much more dramatic 
than the exponential law, which is mainly based on the simple Markovian diffusion approximations. Preliminary MC 
- results [13 can be fitted much better by the form I? ~ [1 — const(A/5'')Af"]. Thus suggests that a critical disorder, 
Ac ~ 6'^A^~", exists and that for strong disorder A > Ac the chain center of mass is localized under the conditions of 
the simulation. 

Such a dramatic reduction of the diffusion coefficient shows that because the dynamics of the chain becomes indeed 
extremely slow the equivalence of the quenched and annealed problems in the infinite volume limit 0, 0, liol lll| has 
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to be questioned for strong disorder, i.e., A > Ac. As soon as A > Ac the chain center of mass is practically localized 
and the chain may only sample intermediate disorder environment. We should implement in this case the quenched 
averaging, i.e. the chain free energy (or any observable physical quantity) must be calculated first for a particular 
configuration of disorder, only after this the average over the disorder can be taken. In any case, as proposed by the 
ref . (3|, ^1 , the system volume indicated by the factor InQ represents then the actually explored volume. For 

completeness we mention that the MC - simulations ^3 ^-l^o reveal the Rouse modes freezing by showing a plateau 
in the time - dependent correlation function and has been predicted first in the ref . [15| . Motivated by our dynamic 
studies we come here back to the static replica approach and reformulate the theory in such a way, that we are able 
to treat the freezing of the chain on internal scales, which correspond to (Rouse) modes. 

The idea is then to use the indication of the center of mass freezing and propose on these grounds a static theory 
for pinned chains in random medium. Clearly then, the theory has to be formulated in terms of the replica theory 
in such a way that the behavior internal degrees of freedom can be treated within the same framework. Obviously 
more than a formulation which is only sensitive to the overall size of the polymer must then be employed. A more 
refined theory needs then more variables than the overall size of the chain, since the behavior of individual segments 
of the chain needs to be studied in view of their freezing behavior. Nevertheless we can expect that within the static 
approach the freezing and nonergodic regimes show up through the replica symmetry breaking (RSB) scheme ^18..19|. 

We reformulate in this paper the replica freezing theory and discuss the different regimes for a Gaussian chain 
by making use the replica approach embedded in the variational formalism 0, l2ll 123^ . Indeed we are going to 
show that the corresponding free energy can be explicitly calculated in the framework of so - called one-step replica 
symmetry breaking (1-RSB) scenario j3i ^M- Moreover it will be formulated more generally, that annealed and 
the replica symmetric theory can be treated as special cases in a more general framework. Within this scenario all 
replicas are grouped into clusters so that the resulting pattern has only two levels: the intra-cluster overlap has a 
finite value whereas the inter-cluster overlap is equal to zero. The first special case, when the cluster contains only one 
self-overlapping replica, corresponds to the replica symmetrical (RS) scenario and in principle leads back to results 
obtained in ref.0- We will show that for such situations a critical value for the disorder exists Ac ~ iV"^+''/-^, which 
corresponds to the Harris criterion (see ^-nd references therein). For stronger disorder the chain is captured by the 
disorder potential and the chain gyration radius (resulting from the interplay between the cntropical and energetical 
terms) is scaled as Rg ^ (1/A)^/'^^^'*-'. This is basically a background for the dynamic localization of the center of 
mass of the chain and the chain collapse to its localization radius Rg. 

In the general case of 1-RSB scenario the free energy is a functional of the mean - square beads (inter-replica) 
deviation within a frozen state V^. This corresponds to an order parameter which characterizes the freezing of the 
internal degrees of freedom inside the chain. We will show that the solution for D"^ arises through the 1-st order 
phase transition (first as a metastable state) at some critical disorder Aa which is similar to critical temperature Ta 
discussed in ref. |2ll 123 . l2 ^ . It is of interest that for reasonably long chains Ac < A^, so that the center of mass is 
localized first and then with increasing A internal degrees of freedom become frozen. 



II. REPLICAS AND VARIATIONAL APPROACH 



A. Model 



We consider a Gaussian chain which is characterized by the d-dimcnsional vector - function R(s), where 1 < s < N 
and labels beads of chain. The chain is embedded in a quenched random potential ,y{R(s)}, so that the whole 
Hamiltonian take the form 



N 



N 



(2.1) 



s=l 



where A = d/2b^ , b is the Kuhn segment length, N is the length of chain and the finite difference VsR(s) = 
R(s -I- 1) — R(s). We will use below the units of measurement where ksT = 1. The quenched random potential 
l^{R(s)} is assumed to be Gaussian distributed with a short-ranged correlator 



(y(r)F(r')) = 



A 



(27ra2)'i/2 



exp 



(r- 



.^2 



2a2 



(2.2) 



where the dispersion A is the main control parameter of the problem. The correlation length a in ea. H2.2fl plays the 
role of the spatial resolution or a minimal scale length. The second line in ea. (|2.2|) corresponds to the limit a 0. 
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B. Replicated variational method 



As pointed out in earlier for the quenched disorder problem the free energy must be averaged over the random field 
F{R(s)}, i.e. [F]^^ = — [InZ]^^, where [. . .]^^ stands for the average over the quenched field and Z is the partition 
function. 

In order to average In Z over V we use the replica trick |l8l Il9| . After introducing n copies of the chain and 
averaging over V with the correlator H2.2|l we obtain 

iFl^^^Um^^^^^^Mm^ , (2.3) 
where the replicated partition function reads 

N . n n N N 



a=l 



with the effective Hamiltonian having the form 



(2.4) 



/ ■n f 71 i\ . 'ri Ti i\ i\ 

n D^s) cxp -ay^Y. [^s^a{s)f + y E E E E (^^'^(■'^) - 

a=l ^ a=l s=l a=l b=l s=l s' = 1 

and Ra(s) with a — 1, 2, ... 7i is the replicated vector - function. The effective free energy 

^ n 

Feff = - In [Z"]^^ = - n ^^-(«) {-^cff} (2.5) 



n N . n AT 

a— 1 s—1 a,fc— 1 s,s' — l 



Because eq. (|2.4ll is not amenable for the exact calculation we use the replicated variational method |2ll l22l |23| . This 
method is based on the inequality 

^^var = i^tr + (i^eff - ^tr) > i^cff , (2.7) 

We take the trial Hamiltonian H^^. in the following bilinear form 

n N n N N n 

^tr = A E E [V.Ra(s)]' + i? E E r^'(^) + CE E i^-(') - (2.8) 

a— 1 s—1 a— 1 s — 1 s — 1 a^b 

with the trial free energy 

/n 
l[D-Ra{s)exp{-Htr} (2.9) 
a=l 

It is worthwile to note that because of the presence of the second term in ea. (|2.8() the "trial" Hamiltonian does not 
preserve the translational symmetry, Ra(s) Ra(s) + C, which is respected by the effective Hamiltonian (|2.6|l . 
Nevertheless it does not damage our results because the chain ultimately is captured by the disorder potential at 
A > Ac in some part of the space (see Sec. HI A). Since this part of the space is not specified one should average 
over the chain locations or (which is equivalent) over the random field realizations (see e.g. Sec. 2.1 in ref.[23|). We 
might expect that the optimization of i^var with respect to B and C leads to a good approximation for Feff, i-e. the 
optimized F*^^. « F^s- After that with the use of ea. H2.3|) we can estimate the free energy [F]^^ = lim„^o F*^^/n. In 
ea. (|2.8|) i? is a variational parameter conjugated to the chain mean-square gyration radius i?g which is responsible 
for the collapse transition. The second variational parameter, C, is conjugated to the mean - square inter-replica 
deviations which characterizes the freezing transition. The form of the Qab - matrix is related with the pattern 
of 1-RSB scenario which we have assumed here. This pattern looks as follows: all n replicas are divided into n/m 
clusters each of which has the size m, so that qab is 1 if a and h {a ^ h) belong to the same cluster and otherwise. 
The averaging < • • • > in ea. H2.7|l means the expectation value with the trial Hamiltonian i?tr- 
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C. Free energy calculation 



The variational free energy i^var can be represented in the form 

Fvai — Fq + Fl ; 

where the entropy term _Fo is 

/ n N 



(2.10) 



(2.11) 



and the interaction part Fi reads 



N 



,.b—l s,s' — l 



(2.12) 



Let us calculate first -Ftr- For this purpose we must first diagonalize the quadratic form (I2.8|l with respect to the 
replica index. The coefRcient matrix of the last part of the quadratic form H2.8|l has the to x to - blocks over diagonal. 
It can be seen easily that this matrix has (n/m) - modes with the eigenvalue and {n — n/m) modes with the 
eigenvalue 2to |2^ |23, 123 ■ Then the eq. I|2.9|l take the form of Gaussian integral 



/n ( n N "S 

n DHais) exp <^ - 5] E [V.Ra(s)]' + Bnl{s) + CkaRl{s)\ \ 
a=l I a=l s=l J 



(2.13) 



where Aa = for 1 < a < {n/m) and = 2to for (n/m) + 1 < a < n. By making use the results of calculation for 
the harmonic oscillator functional integral [23| we have 



Fu 



2 \m) 



In 



sinh(A+) 



sinh(7VA+) 



(---) 

V TO/ 



In 



sinh(A_) 



sinh(7VA_) 



{N - l)dr. 



In A 



(2.14) 



where A+ and A_ can be found from relations 



1 B + 2mC 



A 



= sinh 



(2.15) 



Calculation of (i?tr) can be made by using the scaling arguments |2^. Because Ra(s) are dummy variables -Fir is not 
changed if we scale Ra(s) ^ ^/z'R'^{s) and (i?tr) = [{9/dz) Ftr]^^i. But according to ca. H2.15|l A+ and A_ does not 
depend from the rescaling of coefficients A, B and C. The only part of Fti-{z) which depends from z has the form 
(iV - l)dn\n{zA)/2. After that we have 



{H,,)=n-{N-1) 



For the calculation of 



N 



A{J2Y.lVMs)r)=A—F,, 



(2.16) 



(2.17) 



we make use ea. (|2.14f) . The results reads 

I n N \ 



^ ( EE [^^I^'^(^)]') -"9(^-1) - (-) 2^^Mh[A+]-iVcoth[7VA+]) 



\a— 1 s— 1 



d .dX- 



mJ 2 dA 



(coth[A_] -iVcoth[A^A_]) 



(2.18) 
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For the very long chain, N > 1/A±, sinh(iVA±) « (1/2) exp(A^A±), coth(iVA±) « 1 and eqs. ifnijl - (|??TH|) can be 
substantiaUy siniphficd. For example the trial free energy (|2.14|) takes the form 



dN 



Ft,= -]—X++ n--]—\ 



mj 2 



dN . 



mJ 2 



After straightforward calculation the entropy term H2.11|l becomes 



^0 = (-) 



dN 



mJ 2 



iA 



+ n[l 



1 \ dN 



ml 2 



A_ 



4A 



B+2mC 



(2.19) 



(2.20) 



The calculation of the interaction term Fi (see is more complicated because it involves the replicated Green 

function of the bilinear form l|2.8|l . Employing the exponential representation of the S - function in l|2.12(l we have 



where the replicated correlator 

Qabis, s'; k) 



A 
2" 



E E /7|^^Q-(^'^';k) 

a,6=l s,s' = l ' 



(2.21) 



exp{^k [Ra(s)-R6(s')]} 



exp 



[Gaa{s,s)+Gbb{s',s')~2Gat{s,s')] 



In ea. (|2.22|l the Green function 

Ga6(s,s')= [W-1],JS,S') 

is the result of inversion of the coefficient matrix in the quadratic form 



N 



(2.22) 



(2.23) 



(2.24) 



■d—l s,s' — l 



The explicit calculation of Qab and Fi is given in the Appendix A. It turns out, that the Qab take different forms in 
the respective elements of the matrix. The diagonal elements Qaa read 



Qaa{s,s';k) = expj--^ [5(s,s) +5r(s',s') - 25r(s,s')]| 



exp 



1 - e 



fc^(TO — 1) 

4Amsinh(A+) " J 4Amsinh(A_) 

The off - diagonal elements of Qab inside the blocks are calculated to be 

Qabis, s';k) - cxp|-|^[g(s,s) + g(s',s')-2g(s,s')]| 



1 - e"^-l"""'l 



(2.25) 



exp 



fc2 



1-e-^+l 



4Ar7i sinh(A+ ) L J 471™ sinh(A_ 

Finally the off - diagonal elements of Qab outside the blocks can be written as 

U2 



m — 1 + e 



(2.26) 



Q>(s,.';k) = exp|-|^[g(,s,.)+g(.',.')] 



1 



m — 1 



msinh(A+) msinh(A_) 



(2.27) 



In eqs. (|2.25|l - (|2.27|) g and g stand respectively for the diagonal and off - diagonal elements of the Green - function. 
In addition we have to see how many elements Qab of each type exist for the present problem. It is easily seen that 
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the number of diagonal elements is n. The number of elements is n(m — 1) and similarly the number of elements 
of the Q^f, is n{n — m). 

The final form of Fi is given by 



N 

E 



— A _ \s — s 



, -1 -d/2 



msinh(A+) m sinh(A_) 



+ n{m — 1) 



+ n{n — m) 



1 - e--^+ 



msinh(A^ 



- 1 + 6-^-1''-'' I 
msinh(A_) 



, -1 -d/2 



1 



m — 1 1 



-d/2 



msinh(A+) m sinh(A_) 
and the complete variational free energy is then given by the sum of all the individual contributions. 



(2.28) 



D. Order parameters 



As have been mentioned in the Introduction there are at least two types of phenomena which we face while the 
disorder strength A is growing. At some critical Ac the chain is captured by the disorder potential, so that the chain 
center of mass becomes localized and the mean - squared gyration radius i?g is defined only by disorder. In this case 

i?g play the role of the order parameter and can be expressed in terms of variational parameters B,C as follow 



/UN 
\ c—1 s,s' — l 




A - 1 d 
I ~ nN dB 

TO — 1 



1 



{B + 2mC)Jl 



4A 



(2.29) 



B+2mC 



where we have used eas. (|2.15|l and H2.19|) . In the Sec. Ill we will show that this quenched problem can be treated 
within RS - scenario, i.e. at m = 1. 

On the other hand within 1-RSB scenario the (meta)stable solution for the variational parameter C appears at 
some critical A a- Parameter C is conjugated to the mean - square beads deviation within a frozen state which 
can be written as 



.6=1 s=l 



1 



d 



n{m- l)N dC 



{B + 2mC)Jl 



4A 



(2.30) 



B+2mC 



The order parameter is an inter-replica value which characterizes the mean - squared beads localization 
length and freezing. On the other hand, can be also interpreted as the long - time limit correlation, i.e. 

P2 = limt^oo N-^ ([^(s, t) - R(s, 0)]^). We will consider this in Sec. IV. 



III. CENTER OF MASS LOCALIZATION AND COLLAPSE 
A. RS - scenario 



We will now investigate the quenched problem when the chain is captured by the disorder potential and collapse, 
so that the chain size is determinated only by the disorder strength A. The schematic pictorial representation of this 
situation is given in Fig.l. 
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FIG. 1: a) At the critical disorder Ac chain is captured by the random potential in some part of the space, i.e. becomes 
localized, b) The characteristic scales of the problem: a is the correlation length of the random potential (or the spatial 
resolution), h is the Kuhn segment length, 7?g is the chain gyration radius. Scales satisfy the following relations: a <^ b <^ Rg 



For the purpose of calculations let us take the RS - limit (i.e. n ^ and m = 1) in the interaction l|2.28|l and 
entropy H2.20|l terms. Using ea. H2.28|l and switching from the summation to the integration over the contour variable 
s leads us (after using eq.|^Jl) to the following RS - interaction free energy part 



(1) 



RS 



27r''/2 



N ds 



d/2 



d/2 



(3.1) 



where the mean - square distance 



sinh(A+) 



(3.2) 



Note that this form of the polymer radius corresponds to a simple solution of a Gaussian chain in a harmonic potential. 
The chain statistics is still Gaussian on the interval < s < 1/ A+ while the collapsed configurations can be revealed 
on the larger contour length 1/A-|_ < s < TV. As a result (and because of ^ 1) the chain size is completely 

independent of its total contour length N. Note that the eigenvalue A+ in the limit B A (see below) is small and 
reads 

(3.3) 



The free energy (|3.1|l is typical for the variational method and corresponds to the one - loop approximation. For the 
problem under discussion it was obtained first in Q . Unfortunately the integral in eq. I|3.1|l has the divergency at small 
s and must be regularized. The reason for this divergency lies in the contradiction between the Gaussian statistics 
on the interval < s < 1/A+ and the effective segment - segment attraction which shows up on all scales within the 
variational method (see eg. (13.11) ). The integral regularization can be obtained by the reasonable coarse - graining. 
Namely, we should apply our variational method to the string of statistically independent segments of length I 1/A+ 
(see Appendix B for more details). On this coarse - grained level of description the variational method with only one 
parameter B (as opposed to the case when B{s) is a trial s - dependent function is getting self - consistent and 
the integral in ea. (|3.1|) can be regularized. It should be mentioned another approach regarding similar point. 
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RS - entropy contribution , Pj^g , resulting from ea. H2.20|l in the limit m = l,n ^ and taking into account eqs. 
(|2.15|) and H3.3|) . The final expression for the full RS - free energy F^is can be put in a following form 

^RS - ^RS + ^HS 

dN [b an fB 



fe used also eas. (|B3|) - (|B4|) for 

(j^^ we get 



Id , (3.4) 
where we have used also eas. (|B3|) - (|B4|) for Fj^g' and 1^ respectively. For the mean - square gyration radius from 



Minimization of the free energy H3.4|l with respect to B leads to the solution 

'A 



B*c.A[^] . (3.6) 



The resulting becomes 



Rl « 6^ ) (3.7) 



at d < 4. This result has been obtained first in ref.Q. We can meet the condition B A (see ea. ()3.3|) ') if A <C fe'', 
i.e. R'i ^ (see Fig. lb). The free energy at the minimum B = B* reads 



A 



4-d 



Fks{B*} ^ -N i^^ j <0 , (3.8) 

which corresponds indeed to the form the the free energy confinement Fconf oc Nb^/R'^ for localized chains. As a 
result we conclude that that the chain is captured by the disorder potential at some critical disorder provided that 

N{AJb''f^'~'^^l,i.e. 

A, ~ b'^N-'^+^ . (3.9) 

This criterion has been recently obtained by us within the dynamical consideration in ref.'l5l in various representations. 
It corresponds to the usual Harris criterion for chains in disorder. The quenched RS - free energy H3.4|l and the chain 
confinement is shown in Fig. 2. 

The RS - scenario reveals two facets of the problem: the chain center of mass localization (or confinement) which is 
described by the criterion 1)3. 9|l and the collapse with the characteristic chain size given by eq. 1)3. 7|l . Here it should 
be recorded that as distinct from the results of ref. 0, IE ^| the expression for the chain size does not contain the 
factor InQ (where Q is the volume of the whole random medium) since the chain is effectively pinned for A > Ac. 

B. Annealed problem 

In this case, as we have discussed in the Introduction, the chain is fairly mobile and can span the whole system 
during the time of an experiment. The chain experiences all possible quenched field realizations and the replica trick 
is actually not necessary. The corresponding interaction part of free energy results from eq. (|2.28|) at rt = m = 1 . The 
entropy term follows from ea. (|2.2U|l at n = m = 1 as well. The resulting annealed free energy reads [10| 

dN ^ A ^ 1 



sinh(A+)]'^/^ E 7^ . ,. ,nd/2 • (3-10) 



A 2^mb^' ^ ^^^^ [l_e-A,|.-.'|]'^/2 

By making use eas. (13.2(1 and (|3.3() we can substitute the result of double summation by N"^ / Rj^ whereas the entropy 
term leads to the form b'^N/Rj^. As a result we have 
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FIG. 2: RS - free energy as a function of the variational parameter B and 
different disorder strength A/6'*. The confinement is guaranteed by the con- 
dition that the minimum depth is large enough, i.e. N{A/b'^f^^'^-''^ > 1. The 
upper curve demonstrates the entropy term 



This expression was obtained first in ref.0,|3|. The minimization of ea. H3.11|l leads to the following result. At d > 2 
^ann —oo and Rn 0, or by cutting - off we can state that the chain is collapsed to a size of the order of the 

Kuhn length b. At d < 2 the free energy has a minimum at Rn — b (^AN/b'^Y^'''^ ^ 6 in the long chain limit if the 
cut - off is imposed again. The chain crushing down is of course a result of a Gaussian model. It is well known 
that the more realistic chain with the second v > and the third w virial coefficients chain is collapsed at A > w to 
the globule state with a finite density. 

IV. FREEZING OF INTERNAL DEGREES OF FREEDOM 

As shown in Sec. Ill the Gaussian chain is getting immobile or localized at some critical value of disorder Ac. This 
can be treated by making use the RS - scenario. As a relevant question, we can ask what happened with internal 
motions of the chain beads as the disorder A increases further? It is clear that the chain still experiences the free 
energy landscape and can stay in the different minima or the conformational states. The measure of the distance 
between these states is given by the order parameter which discussed in Sec. II D. The parameter 2?^ is used 
extensively for the heteropolymers, where the landscape dominated statistics can be clearly seen |2^. There is a 
natural mapping in the replica approach 0,113 between sets of states and replicas. The 1-RSB - assumption about 
the clustering in the replica space reflects our view of the states subdivision. The general pattern looks as follows. 
We put together in the same cluster all states that are within a certain distance from each other. The distance 
between different clusters is infinitely large. From the expression for the entropy H2.2()|l it can be seen that the fraction 
of states with a distance 2?^ between them is equal to 1 — 1/m. Fig. 3 illustrates schematically a cluster which is made 
up of three states: a, b and c. In this section we will show that within 1 - RSB - scenario upon disorder growing the 
parameter I?^ can take fairly small value. This means that beads are localized in the globule interior, i.e. the system 
is freezed. Indeed, deep in the frozen regime the change of the chain conformation can only take place on the smallest 
length scales, i.e. ^ a^. 
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FIG. 3: States subdivision in the space of conformations. States of the chain 
within 1-RSB - scenario (marked as a, b, c) form a cluster, so that the distance 
between them is equal O^. To find the distance one traces back the points 
a, b, c along the branches until they merge. The chain can move between these 
states. 



A. RSB - free energy 

The frozen phase is characterized by small fluctuations around individual minima. In this case we assume that 
Debye - Waller factor, 2mC, is large, i.e. 

2TOC>yl>B (4.1) 
and the order parameter given by eq. (|2.^^0|l takes a simple form 

p2~_^L . (4.2) 

2mC ^ ' 

In the opposite limit of small 2mC, i.e. at 

2mC <t:B <t:A (4.3) 

A- ~ A_|_ and it can be seen that the free energy functional Fi{n, m) in ea. H2.28|) and the entropy H2.20|l goes back to 
RS - free energy as it should be. 

In the limit H4.3|l the order parameter V"^ reads 

/ (4.4) 

This is a natural result because in the case when the center of mass is localized but the internal degree of freedom 
are still not frozen (because of small C) the beads localization length is of the order of the gyration radius. 

Let us calculate the RSB - free energy Frsb- For the sake of simplicity we should find first of all the expressions 
for the entropy F^gg (see ea. (|2.20() ') and interaction F-^^-g (see ca. (|2.28|) ') terms in the both limits, H4.1|l and (|4.3|l . 
After that we can find a reasonable interpolating expression, which is suitable for the theoretical investigation. These 
calculations are represented in Appendix C. The overall RSB - free energy which consists of the entropy ljC4|l and the 
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interaction (|C14p terms is given by 



where the dimensionless Debye - Waller factor y = 2mCb^. 



B. Freezing transition 

Now we discuss the transition related with the freezing of the internal degree of freedom and based on the RSB 
free energy expression ea. (|4.5|l . Minimization of ea. (|4.5() with respect of y yields the solution 

A\ ^ m2(m^+i - 1) 



r^i^i) — 77T^ (4-6) 



(m - 1) 



where we have used for A+ the results of RS - solution i.e. A+ ~ {B*/A)^/^ ~ {A/b'^Y^'''^ When minimizing 

eq. H4.5|l we also took into account that y* is close to ?/max, i-e. y* < j/max, where j/max is given by 

52 

?/,nax ^ 2mCme.^b^ - ^ > 1 ' (4-7) 

The stationary condition with respect to m, i.e. idFB.sB / dm)^^^ — 0, leads to the optimal result for m = m* 
defined by 

2 

m* ~ [^" ^""4^ > 1 (4.8) 
rA_] (4-£i)(4+£i) 

Solution y — y* appears first as a metastable one at some A = Aa, which is similar to the critical temperature Ta 
in heteropolymers discussed in ref. 22, 23] . This transition in a metastable state is also of the same nature as in the 
p - spin glass model studied in ref.(27|. 

In Fig. 3 we have plotted the RSB - free energy expression ea. (|4.5|) as a function of Debye - Waller factor y at 
the different disorder strength; in doing so we have used for m the optimal value given by ea. (|4.8() . The metastable 
solution shows up at Aa = 0.753, so that y* grows with the disorder until it hits the ymax, i-e. y* — j/max = b^ /a^- In 
our case it is happened at Ac/b^ ~ 0.796 with t/max = 66.3. At this point all beads are totally localized or frozen. 

It is interesting to make a link with the heteropolymer theory [2^ |2^ . At m* 1 the RSB - free energy -Frsb — 
and eg. 14.81) immediately leads to the corresponding expression for A = Ak which meets this condition, i.e. 

^^[Inj/n.ax]'^ . (4.9) 

This expression quantitatively corresponds to the Kauzmann temperature Tk in heteropolymers |22l |. In ref. |22l| it 
was found that a/ksTK (where a is the variance of the random interaction) is proportional to the square root of the 
entropy loss, which is similar to ea. (|4.9|) . The dependence of y* from A at the different t/max = b'^/a^ is shown in 
Fig. 4. The angular points, where y* hits t/max, correspond to A = Aq where the chain becomes fully frozen. 

In the disorder interval Aa < A < Aq the chain is in a metastable state and the beads localization length lies in 
the range < V'^ — ^g- "^^^ beads become localized when the effective barrier height between the minimum F* 
and maximum F** in Fig. 3 is large enough to assure the beads confinement. Fig. 5 shows the result of the numerical 
estimation of the effective barrier height F^/N = (F** — F*) /N. The result of the fitting can be written as 

i^V^ = 2.03(A- Aa)^-"*^ . (4.10) 
The beads confinement condition, F^ > 1, leads to the following finite chain scaling for the transition disorder 

^A(iV)-AA + ^ (4.11) 

The ea. (|4.11|l means that for the reasonably long chain the internal freezing does not depend from the chain length. 
The comparison with the scaling low eq. 1)3. 9|l for the center of mass localization critical disorder Ac leads us to the 
conclusion that Ac < Aa at A*" ^ 1. By this is meant that as the disorder increases the chain center of mass is localized 
first at A = Ac and then at A = Aa the internal degree of freedom start to freeze. Furthermore at A = Aq > Aa 
the whole system is freezed out with the beads localization length ~ a^. 
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Debye- Waller factor, y 

FIG. 4: RSB - free energy at the different disorder and fixed j/max = 66.3. 

The arrow shows the minimum position y* when it hits ymax, i.e. y* — 2/max. 
The corresponding disorder strength A = Ag = 0.796&^ 



V. CONCLUSION 



We have studied the behavior of the Gaussian chain in a quenched random short correlated field which is specified 
by the field dispersion A. The consideration is based on the variational method in replica space which was previously 
used for the heteropolymer problem [2ll 123 . |2^ within the simplest 1-RSB - scenario. This assumption allows to 
calculate the general free energy functional in an explicit form and consider its RS and RSB cases. 

We have shown that the RS - limit (n — > and m = 1) allows to treat the center of mass localization at some 
critical Ac (see ea. ()3.9|l ^ as well as the chain collapse (see ea. ljB.TIl V In line with ea. (|3.9|l it was found by the methods 
based on the analytical Langevin dynamics calculations 15] that the localization threshold for the center of mass of 
the self - avoiding chain is scaled as Ac ~ }jdj<^-2+vd^ where v is the Flory exponent. More precisely, for the problem 
in question there are two characteristic time scales. One of them is the time of the chain capture tcap by the disorder 
potential. Within our dynamical calculation based on the Hartree approximation |0| (which in turn is reminiscent 
of the mode - coupling approximation T^i) icap — *■ oo. The other time scale is the equilibration time of the internal 
degrees of freedom tcc^■ It should be realized that our results are valid for icap ^ ieq- The main conclusions of ref.[l5| 
have been recently corroborated by the direct MC - simulations |l7j | which among other things shows a pronounced 
decrease of the center of mass diffusion coefficient in the vicinity of A^. The resulting ea. H3.7|l corresponds to the size 
of this trapped chain. Nevertheless the size itself is not sufficient to describe the physical situation. Here we added 
the discussion of the internal degrees of freedom by dividing the replicated system into n/m blocks. 

Then the 1-RSB - limit (i.e. at n ^ and m = m* , where m* is an optimized value of the breaking point parameter 
m) of the free energy expression reveals the nontrivial solution for the inter- replica value V^, the beads localization 
length. This solution appears discontinuously at the critical disorder A = Aa which corresponds to the freezing of 
the chain internal degree of freedom. Beads localization length V"^ falls to ~ with the disorder increasing in 
the interval A a < A < Aq- These results are in a qualitative agreement with the free zing scenario in heteropolymers 
[22', '23]. On the other hand, according the findings based on the dynamical approach Toj the nonergodicity function 
shows up (at least for the long Rouse modes) continuously. Moreover the critical threshold for internal Rouse modes 
freezing is scaled as Ac ~ If'-N"^-'^^ which at first sight can not be reconciled with 60. 1)4.11(1 . At the moment it is not 
clear what is the reason for this discrepancy. Presumably by going beyond the scope of 1 - RSB - scenario to the full 
scale RSB - hierarchy we could fill the gap between dynamical and replica methods. 
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FIG. 5: Dependence of the minimum position y* from the at the dif- 

ferent i/max- The angular point corresponds to the disorder, Aq, when the 
localization length ~ a^, i.e. the chain becomes fully frozen. 



Barrier-fit: 2.03 (A -Ay^)'-*^ 

Barrier: o 



A^ / = 0.753 




0.76 



0.78 



0.8 



0.82 



0.84 



Disorder strength, A / b 



0.86 



FIG. 6: The barrier height as a function of the disorder. The scaling form 
fitting corresponds to the ea. l4.1U|l . 
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APPENDIX A: REPLICATED GREEN FUNCTION FOR THE GAUSSIAN TRIAL CHAIN 



Diagonalization of the quadratic form (|2.8I) in the replica space leads to ea. (|2.13|) . The corresponding Green 
function has two components, G^{s, s'), which associated with two types of eigenvalues. The Green function satisfies 
the equation of motion 



-V,V, + (c^±)']g±(s,s')=^s.' , 
subjected to free ends boundary conditions 



(Al) 



(A2) 



In eqs. (jSI} - (ESj the finite differ ences are defined as follows: VsG^(s, s') = G*(s + 1, s') — G^(s, s'), VsG*(s, s') 
G^(s, s') — G*(s— 1, s'); the "characteristic frequencies" = \J Bj A and lo" 
of eq. HAip becomes the form |^ W\ 



yJ(B + 2mG)/A. The exact solution 



G±(s,s') = 



■ cosh[A±(s-i)] cosh[A±(Ar-s' + i)] 
sinh^Aj-] sinh[7VA±] ' 



if s < s' 



cosh[A± (s'-i)] cosh[A± (iV-s+i)] , 
. sinh[A±] sinh[JVA±] ' U « > 

where A+ and A_ are defined by ea. H2.15() . 

For the large chain, 3> 1/A±, ea. l|A3|l can be written in a much simpler form 



(A3) 



G±(s,s') 



1 



2sinh(Aj 



■exp{-A±|s - s'|} 



(A4) 



In order to calculate the Green function matrix elements Gaf,(s, s') we recall that the n x n - matrix has the block 
- diagonal structure with each block of size m. Let us denote all diagonal elements Gaa =5, off - diagonal elements 
which belong to blocks Gab = g and all elements outside of blocks are equal to zero. For the purpose of diagonalizing 
the matrix Gab we diagonalize first each mx m - block. The corresponding characteristic equation for the eigenvalues 
X is 



g-x _9 
g g-x 



= 



g 9 ■■■ g-x 

It is easy to find that the ea. (|A5|) has one solution 

X+ = 5 - (1 - m)g 

and (m — 1) - solutions 

The determinant of Gab is an invariant of orthogonal transformations and can be represented in the form 

det{Ga6} = (x+)(-) 

= (G+)(-)(G-)("^-) , 

which leads to the identity x+ = G^ and X- ~ G^ and allows to wright down 

G+ = 5- (1 - m)g 
G^ = g-g 

As a result the expressions for non - zero matrix - elements read 

1 



(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



[G^ 



G- 



1 



= -G+ 
m 



1 



1 \G- 

m , 



(AlO) 
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Now we are in a position to calculate correlators Qab{s, s'; k) in ea. (|2.22() . The diagonal element Qaa is 



i(s,s';k) = exp 



AA 



{s,s)+g{s',s')-2g{s,s')] 



exp 



1 - e 



k'^{m — 1) 
AAm sinh(A_) 



1 - e 



-A- \s~s I 



AAm sinh(A_ 

where we have used eqs. (jA10|) and (|A4|) . The off - diagonal elements of Qab inside the blocks are 



(All) 



AA 



{s,s) + g{s',s')-2g{s,s')] 



exp 



fc2 



1 — e^-'^+l''^'*' 



4Arnsinh(A+) L ' J 4Ato sinh(A_) 

Finally the off - diagonal elements of Qab outside the blocks can be written as 

U2 



m — 1 + e" 



-A- |s — s I 



(A12) 



Q>(s,.';k) = exp|-|^[g(,s,.)+g(.',.s')]| 



1 



m — 1 



msinh(A+) msinh(A_) 



(A13) 



How many elements Qab of each type exist? It is easily seen that the number of diagonal elements is n, the number 
of elements is n{m — 1) and finally the number of elements is n(n — m). After that the interaction part of 
free energy (see ea. (|2.21ll ) can be represented in the general form 



Fi{n,m) 



N 

E 



A 



+ n(m — 1) exp < — 



n exp 
fc2 



fc2 



4:Am sinh(A 



+ n(n — m) exp 



„2 r 



AA 



A Am sinh(A 

m — 1 



I _ g-^+k-s'l 



k'^{m — 1) 



4Am sinh(A_) 



AAm sinh(A„) 

TO - 1 + 6"^-!''^"' 



1 - e-^-l''"'*' 



TOsinh(A+) TOsinh(A_) 



(AM) 



or after integration over k it can be written as 

N 



Fi{n, to) 



A 



+ n(m — 1) 
+ n{n — to) 



E 



l_^-\+\s-s'\ i_e- 



A_ Is — s 



-| -d/2 



1 - e 



TOsinh(A+) 771 sinh(A_) 

-1 -d/2 



— A+ Is — s' 



TO sinh(A^ 
1 



^- 1 + e--^-l'''-'' I 
TO sinh(A_) 



TO — 1 1 



TOsinh(A+) TO sinh(A_) 



-d/2 



(A15) 



APPENDIX B: REGULARIZATION BY THE COARSE - GRAINING 

First of all it is pertinent to note that the chain statistics is Gaussian on the interval < s < 1/A+ (see ea. H3.2f) ') 
whereas the collapse shows up on the larger contour variable, 1/A+ < s < N. The whole collapsed chain can 
be seen as made up of the statistically independent segments of length I « 1/A+. The number of such segments 
ni = N/l = iVA+ >• 1. The statistical independence and the equipartition of energy imply that the free energy is 
proportional to fc^Tn/, i.e. 

2 

^^ksTN f^j (Bl) 
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and we go back to eq. H3.8() . Now we can consider the chain as a string of totally collapsed statistically independent 
segments of length / « 1/-^+ • This coarse - grained picture provides a basis for the regularization of the integral in 
eq. (jSH. 

The reason for the divergency of this integral at small s is the following. There is a mismatch between the Gaussian 
chain in the external field B and the variational method where the value of B is expected to find from the variational 
free energy minimization. The divergency on the interval, < s < 1/A+ shows that the effective segment - segment 
interaction within the variational method comes in a contradiction with the Gaussian statistics on this interval. One 
way to resolve this contradiction is to take the trial free energy in a more general form, X^fLi R-(s)i?(s — s')R(s'), 
i.e. now B{s) is a trial function but not a constant field. That actually what we have considered in the recent papers 

Nevertheless, at the moment it is enough to stay with only one variational parameter B which correctly describes 
the behavior of the string of statistically independent segments. In other words, the variational method should be 
applied to the coarse - grained model which we have discussed above. 

Mathematically this simply amounts to the limiting of the contour variable s on the interval 1/A+ < s < oo, i.e. 
eq.(|OJ reads 



5.(1) 
RS 



A 



After that the RS - free energy becomes 



iV[sinh(A+)]''/^ j ds 

1/A+ 



1 



1 - e-^+' 



(B2) 



^(1) ^ _ 



A ^[sinh(A4 



]d/2 



(B3) 



where 



Id = dx 



(B4) 



so that the integral in ea. (|3.1|) is regularized. 



APPENDIX C: RSB - FREE ENERGY FUNCTIONAL: LIMITS AND INTERPOLATION 



Let us calculate first the entropy term (|2.20() in the limits 14. 1() and (|4.3() . For the eigenvalue A_ (see ea. (|2.15|) ') one 
can write 



A_ = 2 In 



1 /S + 2mC 



1 



B + 2mC 



(CI) 



A V iA 

i.e. at C = A_ = A+. We expand eas. HCl|l with respect to small C retaining only the main terms. Then the 

5.(0)- 
RS . 

m-1 dN mC 



corresponding RSB - entropy term ( after subtraction of RS - term, i.e. -F^gg — [hm„^o 

Fo/n]-F^^) reads 



-^RSB 



m 2 AA 

For the large C using again ea.HCl|l in ea. H2.20|l yields the result 

m-ldN 



(C2) 



(0) 



RSB 



" /2mC\ 




rl A ) 


- 1 







(C3) 



The suitable interpolating expression which embraces both, eas. (IC2|l and (|C3p can be written in the simple form 



(0) 



m-ldN 



RSB 



In 



2mC 



A 



(C4) 



In such form the entropy term has been obtained in ref |22l |2 
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Now we can make the same calculations for Fi{n,m) in ea. (|A15|) . Because the dependence from C is contained 
now only in A_ it is convenient for small C to represent 



A+ 1 + 



mC 



(C5) 



Expansion of the first term in ea. l)A15|l has the form 



X _ g-^+|s-s'| m— 1 1 — e^-^- 1*""'! 



m sinh(A_| 



sinh(A_ 



-d/2 



-d/2 



I — g-^+k-s'l 

sinh(A+) 

A+|.s^s'|e-^+l^-"'l 
1 — 6--*'+!'*-'*' 

By the same way the 2nd and 3rd terms in eg. ljAlSjl yields correspondingly 



d m — 1 mC 
2 ~~m W 

A+ coth(A+) 



1 _ Q-^+\s-s'\ 

msinh(A+) 



-d/2 



(C6) 



msinh(A+) 



sinh(A_) 



-d/2 



[sinh(A+)]^/V^[sinh(A+)]'^/^!^ 
A+ coth(A+ 



— <i'|p--^+N-s'l 



A+|s — s'|e 



m — 1 + e 



— |s — s' 



(C7) 



and 



m — 1 



msinh(A+) m sinh(A_) 



[sinh(A+)]'^/2 + ^ [sinh(A+)]''/2 A+ coth(A+) (C8) 

2 TO is 



The 3rd term in ea. (jA15|l has the factor n{n — to) which at n — > takes the form n{n — to) w —n — n{m — 1), 
so that we can combine the 3rd term with the 1st and 2nd terms correspondingly. After that by switching from the 
summation to the integration over s and s' we arrive at the following two contributions to the RSB - interaction free 
energy part 



p(i) _ 



lim 



Fi{n, to) 



Frs = Ji + J2 



(C9) 



where 



Ji 



AdN TO - 1 toC [sinh(A+)] 



d/2 



A. 



A+ coth(A+) 

d 



Id 



AdN TO-1 toC [sinh(A+)] 



d/2 



2TTd/2l)d ^ 



B 



A. 



[A+coth(A+) + 1] , 



(CIO) 



where Id is given by ea. l|B4|l . At small A+ it can be seen that A+ coth(A+) ~ 1 + X\/2 and both terms, Ji and J2 are 
negative. Then the resulting expression for F^J^ at small C is given by 



^(1) _ AdN TO- 1 toC [sinh(A+)]'*/^ 



RSB 



(Cll) 



2nd/2bd m B A^ 

In the limit of large C the eigenvalue A_- is large also and all terms which include 1/ sinh(A_) can be neglected. 

RSB [linin^o Fi{n, m)/n\ - F^^ 



After that from eo. ljAlSp for -F^sb ~ [^^'^'i-n^Q Fi{n, m) /n] — F^^ we immediately have 



^RSB ~ 



2T:d/2^d 



(to^^^-i) 



Let us denote 



2toC 

y = — — = 2mC¥ 



(C12) 



(C13) 
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Then the parameter mC/B ~ y/(2A^), where we have used ~ B/A. After that the expression which interpolates 
between eas. (ICll|) and (IC12II can be written in the form 



^(1) ^ AdN [sinh(A+)]'^/^ imi+^-l) y 

^RSB - 2^,/2;,d m(m^ + ^-l) ^ ' 
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